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Localized quantum wave packets can be produced in a variety of phys- 
ical systems and are the subject of much current research in atomic, 
molecular, chemical, and condensed-matter physics. They are particu- 
larly well suited for studying the classical limit of a quantum-mechanical 
system. The motion of a localized quantum wave packet initially fol- 
lows the corresponding classical motion. However, in most cases the 
quantum wave packet spreads and undergoes a series of collapses and 
revivals. We present a generic treatment of wave-packet evolution, and 
we provide conditions under which various types of revivals occur in 
ideal form. The discussion is at a level appropriate for an advanced un- 
dergraduate or first-year graduate course in quantum mechanics. Ex- 
plicit examples of different types of revival structure are provided, and 
physical applications are discussed. 



Accepted for publication in the American Journal of Physics 
1 Address after September 1, 1995: Physics Department, Cornell University, Ithaca, NY 14853. 



I. INTRODUCTION 



The advent of short-pulsed lasers has made it possible to produce and detect co- 
herent superpositions of quantum-mechanical electron states for a variety of physical 
systems. Once produced, such superpositions can result in the formation of localized 
electron wave packets. The evolution and dynamics of these wave packets are the 
subject of much current investigation in many areas of physics and chemistry. 

One of the most studied applications of wave packets, both theoretically and 
experimentally, is in atomic physics using Rydberg atoms. Excitation by a short laser 
pulse produces a Rydberg wave packet, i.e., a superposition of highly excited single- 
electron states [1], 0. One aim of investigations on Rydberg wave packets is to probe 
the interface between classical and quantum mechanics M. This is feasible because, 
when a Rydberg wave packet first forms, its motion is periodic with the same classical 
period T c \ as a charged particle in a Coulomb field @, |]. However, this motion lasts 
for only a few cycles, whereupon quantum-interference effects cause the wave packet 
first to collapse and then to undergo a sequence of revivals fj], § || [7], ||, ||, [K| |TTJ . 

The revivals of Rydberg wave packets are characterized by the recombination of 
the collapsed wave packet at a time t mv into a form close to the original shape that 
again oscillates with period T c \. This may be viewed in part as a manifestation of 



the quantum recurrence theorem [T^, [L3|], which follows because the time evolution 
is driven by discrete energy eigenstates. In addition to the full revivals, fractional 
revivals occur at times that are rational fractions of the revival time t mv and corre- 
spond to the formation of macroscopically distinct subsidiary wave packets pi. These 
have periodicity equal to a rational fraction of the period T c \. Both full and fractional 
revivals have been detected in pump-probe time-delayed photoionization experiments 
14] , 15 , |T6| , 17f1 . More recently, it has been shown that, for times beyond the revival 



time t rev , a new sequence of full and fractional revivals commences || |K| |ll[]. These 
are characterized by a longer time scale called the superrevival time t sr . At times that 
are certain rational fractions of t sr , distinct subsidiary waves again form, but with a 
periodicity that is a rational fraction of the revival time t rev . These long-term frac- 
tional revivals culminate with the formation of a single localized wave packet, called 
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a superrevival ||, that more closely resembles the initial wave packet than does the 
full revival at time t rev . 

The study of revivals is not limited to Rydberg wave packets in atomic physics. 
Wave-packet dynamics is being explored in several areas in physics and chemistry, 
and the structure of revivals is relevant to a variety of problems. We list here some 
examples. In molecular physics, a new realm of phenomena involving wave packets 
has opened up with the emergence of femtosecond pulse technology ITS]]. Ultrashort 
pulses are used to prepare and probe molecular wave packets in ways that permit the 
charting of chemical processes. A related application attempts to control chemical 
reactions using optimized light pulses |L9|, |20]]. For example, frequency-chirped laser 
pulses have been synthesized to control the evolution of vibrational wave packets of 
iodine Wave packets have also been produced in semiconductor quantum- well 



systems [22]. Their dynamics is similar to that of atomic systems, allowing for the 



observation of quantum beats and revivals [[23|, Revival structure has also been 
observed experimentally in the evolution of the atomic inversion in a micromaser 



25|j . Some recent theoretical predictions include fractional revivals for the Jaynes- 
Cummings model ||26| , for wave packets in a Morse-like potential |2"7 ], and for an atom 
bouncing vertically in a gravitational cavity 



Several different types of revival structure can be found in these systems. Most 
of them share the feature of exhibiting an initial periodicity, often corresponding to 
that of a related classical problem, followed by a sequence of collapses and subse- 
quent revivals. However, the times at which the revivals occur, the extent to which 
they resemble the initial wave packet, and the question of whether fractional revivals 
and/or superrevivals occur all depend on the quantum system being considered. 

In this paper, we show that the revival structure of wave packets for different 
quantum systems can be treated in a generic fashion and that classes of quantum 
systems exist for which the revival structures are perfect. We define time scales 
that govern the behavior and evolution of a quantum wave packet and use them to 
introduce the various types of revival that can occur. The results are illustrated 
with examples of simple quantum systems, including the harmonic oscillator, the 
infinite square well, the rigid rotator, and Rydberg wave packets. In addition to 
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being standard solvable systems in quantum mechanics, these examples also serve 
as the foundation for much of the theoretical analysis performed on actual physical 
systems of interest in atomic, molecular, chemical, and condensed-matter physics. 
Other examples could be developed via a straightforward extension of the methods 
presented here. 

The level of treatment and the examples we consider make this topic appropriate 
for an advanced undergraduate or first-year graduate course in quantum mechanics. 
The superposition of eigenstates that comprises a wave packet and the analysis of its 
time evolution can be studied either using standard software packages for symbolic 
manipulation and graphics or via computer programs comparable in difficulty to 



others successfully employed by undergraduates [f2^, . The treatment of Rydberg 
atoms presented here is particularly useful in discussions of the classical-quantum 
correspondence and in much current research in laser-based atomic physics. 

This paper is organized as follows. In the next section, we describe the generic 
features of localized quantum wave packets, define the time scales that control their 
evolution, and discuss the revival structure. Specific examples are considered in Sec. 
III. The autocorrelation of a time-evolved wave function with the initial wave function 
provides a striking illustration of different types of revival structure. This is described 
in Sec. IV, along with illustrative examples. We summarize in Sec. V. Throughout this 
paper, we work mostly in atomic units for which h = e = m e = 1. The conversions 
to standard units are 1 a.u. ~ 5.29 x 1CT 11 m for lengths and 1 a.u. ~ 2.42 x 1CT 17 
sec for times. 

II. REVIVAL STRUCTURE OF QUANTUM WAVE PACKETS 

A. Generic Features 

The time-dependent wave function for a localized quantum wave packet formed 
as a one- dimensional superposition of energy eigenstates may be written as 

*(f, t) = J2 Cn^„(r) exp \-iE n t] . (1) 

n 

Here, the vector r denotes the position of the wave function in one or more dimensions. 
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The label n represents the quantum number relevant for the sum in the superposition. 
The wave functions if) n (r) represent (possibly degenerate) energy eigenstates with 
corresponding eigenvalues E n . Any quantum numbers other than n are suppressed 
for convenience. The coefficients c n are the weighting factors in the superposition, 
given in terms of the initial wave function by c n = (ip n (f)\^?(f, 0)). 

In the general case the superposition (fl|) may include continuum states, for which 
the index n would become continuous and the sum would become an integral. How- 
ever, our primary interest here is in superpositions of bound states with negligible 
continuum contributions, as these are the type produced in typical experiments on 
revivals in wave packets. Therefore, in what follows we assume the index n is discrete. 
Also, we assume that the expansion ([!]) is strongly weighted around a mean value n 
for the quantum number n. This assumption is reasonable for situations involving 
the production of a localized wave packet using a short-pulsed laser. The uncertainty 
principle applied to the laser light states that tAE > 1/2, where r is the duration 
of the laser pulse and AE is the spread in energy. Exciting a quantum system from 
its ground state with a short laser pulse produces a superposition of states centered 
on a value n that depends on the mean frequency of the laser light and has a spread 
a in n that is inversely related to the duration r of the laser pulse. In situations of 
interest n is typically large, so we assume this in what follows. 

We model the weighting probablities |c„| 2 as a gaussian distribution 

I (n-n) 2 

e 5^ . (2) 



This simple assumption provides a symmetrical distribution in n with mean n and 
standard deviation a. It therefore can be used to match the mean value and spread 
of states produced by a laser, though the shapes of the distributions may differ. 
The assumption of a gaussian distribution is unimportant in what follows, serving 
primarily to simplify computation of the wave-packet evolution. Other weighting 
distributions could also be considered. For example, coherent states for the simple 



harmonic oscillator [31] may be produced by weighting the energy eigenstates with 
the coefficients 

c n = e-^ 2 ^= , (3) 
vn! 
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which depend on a complex parameter a. This distribution is asymmetric in n and is 
weighted around n ~ | ct | 2 with spread a ~ |a|. These coherent states have minimum 
uncertainty, exhibit motion with the classical period, and maintain their shape with 
time. Minimum-uncertainty coherent states may also be obtained for other problems 
32[ , although the associated wave packets subsequently collapse and undergo revivals. 

The wave packets produced in experiments using single short laser pulses to excite 
Rydberg atoms are localized only in the radial coordinate and have p-state angular 
distributions. They are called radial Rydberg wave packets. Minimum-uncertainty 
solutions have been obtained for these wave packets |34|. They are examples of 
squeezed states, which are characterized by having a shape and uncertainty product 
that oscillate with time []35|] . To produce Rydberg wave packets that move on an 
elliptical trajectory in a plane, external fields must be applied to mix the angular 
states. One scheme for producing elliptical wave packets has been proposed |36|. A 



squeezed-state description of planar elliptical wave packets has recently been obtained 
f37| . The expressions for the corresponding weighting coefficients for the radial and 
elliptical squeezed states are given in refs. |34|, |37 . 



B. Time Scales 

The assumption that the weighting probabilities | strongly centered around 

a mean value n means that only those states with energies E n near the value E n enter 
appreciably into the sum in Eq. ([I]). This permits an expansion of the energy in a 
Taylor series in n around the centrally excited value n: 

E n ^E n + E' n (n-n) + ^(n-n) 2 + ^'(n-nf + --- , (4) 

where each prime on E n denotes a derivative. 

The derivative terms in Eq. (0J) define distinct time scales that depend on n: 

2tt 2tt 2tt 

~~ "TfM ' * rev ~~ i I F"\ ' ^ sr ~~ il F"'\ ' ^ 

I n I 2 I n I 6 ' n ' 

The first time scale, T cl , is called the classical period. The second time scale, t rev , is 
the revival time. The third time scale, t sr , is the superrevival time. In practice, for 
all cases of interest the expansion (HI) is such that T cl t rcv t ST . 
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Keeping terms through third order and disregarding an overall time-dependent 
phase, we may rewrite Eq. ([[]) as 



*(T> t) = H c n ip n (r) exp — 2ttz 



n 



n — n)t (n — n) 2 t (n — n) 3 t 




) 



(6) 



Since both n and n are assumed to be integers, the powers of (n — n) are all integer- 
valued as well. 

A more general treatment of the energies and laser excitation would allow for 
noninteger values n* of the quantum number as well as a laser detuning with E n 
off resonance. The consequences of these additional effects have been considered for 
excitations of Rydberg wave packets in alkali-metal atoms, which have energies E n * 
given by a Rydberg series depending on quantum defects [[|. The quantum defects 
shift the values of n* to noninteger values and thereby model the effects of the core 
electrons in an alkali-metal atom. Ref. shows that the effects of quantum defects on 
the long-term evolution and revival structure of a Rydberg wave packet are different 
from the effects caused by a laser detuning. The noninteger values for n* and the 
laser detuning cause shifts in the revival time scales and periodicities that can be 
distinguished from each other. The treatment of the wave-packet evolution including 
these effects is similar to that given here. Since the time scales are much simpler 
for the case of integer-valued quantum numbers and with the laser on resonance, we 
restrict our analysis here to these cases. 

C. Revival Structure 

The expansion (|6]) shows that the time dependence of the wave function is governed 
by the three time scales in Eq. (|5|), which in turn are controlled by the dependence 
of the energy on the quantum number n. In what follows, we first summarize the 
behavior of a wave packet for the generic case where all three time scales are finite. 
We then discuss special classes of quantum systems for which wave-packet evolution 
is particularly simple. 

For small values of t, the first term in the phase in Eq. (||) dominates. Evidently, 
during this interval motion of the wave function \&(f, t) is approximately periodic in 
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time with period T c \. 

As t increases and becomes appreciable compared to t Tev , the second term in the 
phase modulates this behavior, causing the wave packet to spread and collapse. How- 
ever, at later times near t rev , the second term in the phase in Eq. is approximately 
equal to 2iri, which makes it irrelevant, and once again the motion is governed by 
the first term. As a result, the wave packet regains its initial shape, and its motion 
is periodic with the period T c \. This is called a full revival. For times greater than 
t rev but much smaller than t SI , the revival cycle repeats. At every multiple of t rev , the 
second term in the phase again equals 2iri and a full revival occurs. 

At special times that are rational fractions of t rev , the wave packet gathers into 
a series of subsidiary waves called fractional revivals. The motion of these fractional 
revivals is periodic || with period given by a rational fraction of T c \. 

The presence of the third-order term in the time- dependent phase modulates the 
whole cycle of periodicity, collapse, and fractional /full revivals |TT[. At times appre- 
ciable compared to the superrevival time scale t BT , the full- revival structure collapses 
and a new system of revivals appears with features controlled by t sr . At times given 
by the rational fractions -t ST , with q an integer multiple of 3, the wave packet again 
separates into a series of subsidiary waves called fractional superrevivals. In this 
case, however, the motion of the wave packet is periodic with a period equal to ^rev- 
Moreover, at the time |t sr , the wave function reforms into a single wave packet that 
resembles the initial one better than does the full revival at time t Tey . This new 
structure is called a superrevival. 

We have kept only the first three derivative terms in the energy expansion (|) 
and hence have only three time scales in Eq. (||). In principle, higher-order terms 
contribute and lead to a modulation of the full/fractional superrevivals for times t > 
£ sr . In most physical systems, these higher-order effects are irrelevant. For example, 
for Rydberg wave packets the fourth-order time scale is comparable to the lifetime of 
the Rydberg states and hence the Rydberg wave packets spontaneously decay before 
this time scale comes into play. The superrevival time scale t BT , however, is several 
orders of magnitude less than the spontaneous lifetime as well as the lifetime for 
blackbody-induced transitions, so it is feasible for experiments to detect superrevivals 
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of Rydberg wave packets. 

We can use the above analysis to identify classes of quantum systems that clearly 
exhibit each of the above behaviors. The simplest nontrivial class consists of systems 
which are perfectly periodic and for which wave packets undergo neither revivals nor 
superrevivals. By perfect periodicity, we mean that the motion of the wave packets 
between times t = and t = T c \ exacly matches the motion between times t = mT c \ 
and t = (m + 1)T C \ for any positive integer m. For such systems, there are no revivals 
or superrevivals in the sense that the initial periodic motion is maintained indefinitely. 
In this case, the time scales t sr and t rev are infinite, or equivalently d k E n /dn k \n = for 
all k > 2. This means that wave packets in quantum systems with energy eigenstates 
of the form E n = A + Bn, with A and B independent of n and B ^ 0, exhibit perfect 
periodic motion and regain their initial shape at times that are multiples of T c \. The 
wave packets never fully collapse and there are no revivals or superrevivals. Examples 
of systems of this type are provided by the D-dimensional harmonic oscillators, for 
which the eigenenergies are E n = (n + D/2)u, where uj is the frequency and n is 
the principal quantum number. Since these systems have no higher-order modulating 
terms in the time-dependent phase of Eq. (|5]), the shape of the wave packet remains 
unchanged after each classical period. 

The next simplest class consists of systems for which wave packets undergo revivals 
but not superrevivals. This means t sv is infinite, or equivalently d h E n /dn k \n = for 
all k > 3. In this class are quantum systems with eigenenergies having n dependence 
of the form E n = A + Bn + Cn 2 , where A, B and C are independent of n and C 
is nonzero. The associated wave packets undergo motion with the classical period 
modulated by the revival phase, as expected. However, in these systems the full 
revival is a perfect replica of the original packet, since the revival structure is not 
modulated by a superrevival timescale. The infinite square well and the rigid rotator 
are examples of quantum systems having energy spectra in this class. 

One might suspect the existence of a class of quantum systems for which the 
superrevivals are perfect in the sense that there are no higher-order contributions 
to the time-dependent phase of Eq. (^). This would be the case if there existed a 
quantum system with bound-state energies having a cubic dependence on n. However, 
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the bound-state energies in nonrelativistic systems with confining potentials can rise 
no faster than the power n 2 in the limit of large n |38| . In principle, a quantum system 
might be identified with energies varying as the cube of n for a restricted range of 
n. If a wave packet in such a system were constructed solely as a superposition of 
energies in this range, it would have perfect superrevivals. We are unaware of any 
physical quantum system with this feature. In the remainder of this paper where 
superrevivals are considered, we take the hydrogen atom with E n = —l/2n 2 as our 
canonical example of a quantum system with eigenenergies non-quadratic in n. Wave 
packets in this system exhibit the full range of behavior described above. 

III. EXAMPLES 

In this section, we illustrate the discussion in Sec. II with examples of wave packets 
for several different quantum systems: the harmonic oscillator, the infinite square well, 
the rigid rotator, and the hydrogen atom. 

A. Simple Harmonic Oscillator 

According to the discussion in Sec. II, wave packets for the simple harmonic os- 
cillator never collapse and have no revivals. We illustrate this by plotting the wave 
packet in Eq. ([[]) using the gaussian distribution for the c n coefficients given in Eq. 
(@) and the standard energy eigenfunctions ip n (x) for the simple harmonic oscillator. 
Since the sum in Eq. (||) is strongly peaked on the value n, we may truncate the sum 
after a finite number of terms and numerically evaluate the result. For convenience, 
we set u = 1 in atomic units. 

Figure 1 shows various stages in an orbital period of a wave packet with n — 15 and 
cr = 1.5. The wave packet is initially located at the right-hand turning point of the 
corresponding classical oscillator. Note that, although the wave packet is localized, 
its shape is not gaussian and so it is not a coherent state. The shape therefore changes 
during the motion. Nonetheless, the packet follows the classical motion with period 
T c \ = 27T / E' n = 2ti. The motion is perfectly periodic, so that each subsidiary figure 
shows the wave packet at the given time up to an integer multiple of T c \. The wave 
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packet does not collapse and there are no revivals or superrevivals. 

B. Infinite Square Well 

The energy eigenstates for a particle in an infinite square-well potential with 
< x < L are ip n (x) = y/Jsin 1 ^, and the eigenenergies are E n = "^jt- For 
simplicity, we consider a box of length L — 1 in atomic units. With these values, 
T c \ = -3- and t rev — - = 2nT d for a superposition centered on the value n. The 
discussion in Sec. II shows that the wave packet is initially periodic with period T c \ 
before subsequently collapsing and exhibiting fractional and full revivals. The full 
revivals are perfect revivals, since the superrevival time t sr is infinite. 

Figure 2 shows the time evolution of a wave packet in an infinite square well. 
We have used Eq. (^) for the c n coefficients with n — 15 and a = 1.5. Figure 2a 
shows the initial wave packet at the left side of the box near the origin. It is highly 
oscillatory because it is receiving an impulse from the infinite potential well. Away 
from the walls, Figure 2b, it is smoother since it moves freely. Figures 2c and 2d 
show the wave packet colliding with the right-hand wall and moving back to the left, 
respectively. At time T c \, the wave packet has completed one cycle and is shown in 
Fig. 2e. Note that the shape has changed. The wave packet is collapsing. At time 
t rev , however, the wave packet regains its initial shape. In fact, the shape of the wave 
packet at any time t exactly matches the shape at time t plus an integer multiple of 

^rev 

C. Rigid Rotator 

A quantum system closely related to the infinite square well is the rigid rotator 
in two dimensions. The hamiltonian is H = -^f, where L z is the angular momentum 
and / is the moment of inertia. For simplicity, we set / = 1 in atomic units. The 
eigenfunctions obey periodic boundary conditions and are given by i() n (<f>) = ~^ em4 \ 
and the energies are E n = n 2 /2. The analysis of revival structure in Sec. II predicts 
the same behavior for the rigid rotator as for the particle in an infinite square well, 
since both problems have quadratic energy dependence on n. 
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Figure 3 shows a wave packet for the rigid rotator with n — 15 and a = 1.5. 
Since n is positive, the wave packet moves counterclockwise (Figures 3a, 3b, 3c). The 
motion is initially periodic, but the shape of the wave packet changes as it moves. 
For times greater than T c \, the wave packet spreads and collapses. Figure 3d shows a 
fractional revival. At t = \t rev , the wave packet has reformed into two macroscopically 
distinct packets moving with period \T C \. Figure 3e shows another fractional revival at 
t = |t rev , consisting of three distinct subsidiary wave packets moving with periodicity 
| T d . At t = |t rev , a single wave packet forms (Figure 3f). This wave packet resembles 
the initial wave packet, but it is a half cycle out of phase with the classical motion. 
The full revival at t — t vcv is a single wave packet that exactly matches the initial one 
in Fig. 3a. 



D. Rydberg Wave Packets 

Depending on the excitation method, several different types of Rydberg wave 
packets can be considered, each with different geometry for the orbital motion. We 
consider a superposition of circular states having I = m = n — 1. These states have 
maximal values for the expectations (L z ) and (L 2 ) and produce a wave packet tightly 
confined to the x-y plane and moving along a circular orbit |J9]. We also restrict our 



analysis here to wave packets for hydrogen, which have energies E n = —l/2n 2 . In 
ref. ||11|| , we examine the revival structure of Rydberg wave packets for alkali- metal 



atoms using an analytical supersymmetry-based quantum-defect model |4(J. We find 
that the quantum defects cause additional shifts in the time scales, but the structure 
of the revivals is similar to that given here for hydrogen. 

We consider a superposition of circular hydrogenic states i^nimif) with I = m = 
n—1 in Eq. (|I|) and a gaussian superposition (□) for the c n coefficients. The time scales 
are given by T c \ = 27m 3 , t rev = ^T c \, and t sr = ^pt rev in atomic units. The analysis 
of revival structure in Sec. II shows that a Rydberg wave packet of this type initially 
follows the classical periodic motion of a particle in a circular orbit. After several 
cycles, the wave packet spreads and collapses. At times that are rational fractions 
of t rev , it exhibits fractional revivals. A full revival occurs at time t rev . Fractional 
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superrevivals occur at times that are rational fractions of £ sr . These culminate at time 
|t sr with the formation of a full superrevival. 

Figure 4 shows the time evolution of a circular wave packet with n = 120 and 
a = 2.5. Converting the time to nanoseconds, we have T c \ ~ 0.26 nsec, 

^rcv — 21.0 

nsec, and t sr ~ 1890 nsec. Figure 4a shows the initial wave packet as a function of 
position in the x-y plane. Throughout the motion, it remains at a constant mean 
radial distance (r) w |n(2n + 1) given by the expectation value of r in the state 
with I = m = n — 1. For n = 120, (r) ~ 14460 a.u. When the wave packet returns 
to its initial position at the end of the first orbit, it has spread (Figures 4b, 4c). It 
subsequently collapses and exhibits fractional revivals. Figure 4d shows the fractional 
revival at ^t rev with two subsidiary wave packets. As expected, the full revival (Figure 
4e) at time i rev does not exactly match the initial wave packet. The full superrevival 
at |t sr is shown in Fig. 4f. This wave packet resembles the initial wave packet better 
than the full revival does at time t rev . 

The example considered here uses a relatively large value of n equal to 120 to illus- 
trate the revival structure. Experiments that have detected full and fractional revivals 
of Rydberg wave packets have involved values in the range n ~ 40 - 90, corresponding 
to revival times in the range t rev ~ 1 - 6 nsec. The apparatus in these experiments 
requires using a delay line on the order of a few nanoseconds between pump and 
probe signals to excite and then ionize the Rydberg atoms. A full superrevival at 
time |t sr can be detected with a similar delay line for wave packets with n ~ 25 - 50. 
For example, with n = 40, the full superrevival occurs after approximately 1.3 nsec, 
making detection possible with current experimental apparatus. 

IV. AUTOCORRELATION 

The autocorrelation function for a time-evolved wave function with an initial wave 
function is given as A(t) = (\I/(r, 0)|\I/(r, t)). Its absolute square gives a measure of 
the overlap between the wave packet at time t = and at a later time t. Substituting 
the expansion in Eq. (HD into |A(t)| 2 yields 

l^)l 2 = |El c «l 2e " ii?n4 2 • (7) 
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Using Eq. (0) for the c n coefficients, we can plot |A(t)| 2 directly for the different 
quantum systems. 

The autocorrelation function gives a striking illustration of the different types of 
revival structures. Numerically, |v4(£)| 2 varies between and 1. When a wave packet 
exactly matches the initial wave packet, |^4(t)| 2 equals 1. If the wave packet is far 
from its initial position, then |^4(t)| 2 — 0. Furthermore, the periodicities in |^4(t)| 2 
reveal the periodic behavior of the wave function. In particular, fractional revivals 
and fractional superrevivals appear as periodic peaks in the autocorrelation function 
with periods that are rational fractions of T c \ and t Tev , respectively. 

Figure 5 shows the absolute square of the autocorrelation function as a function 
of time for four examples. The rigid rotator is excluded because its autocorrelation 
function is similar to that of the particle in an infinite square well. 

Figure 5a shows the autocorrelation for the evolution of a free-particle wave packet 
in one dimension. The wave packet consists of a gaussian-weighted superposition of 
harmonic waves, with momentum-dependent coefficients 



We choose a = 2.5 and po = 10 in atomic units. The wave packet starts at the origin 
and moves away, spreading and gradually collapsing. At t — 0, the autocorrelation is 
equal to 1. For times t > 0, it decreases and asymptotically approaches 0. There is 
no periodicity and no revivals. 

Figure 5b shows the autocorrelation for the harmonic oscillator. The wave packet 
is readily seen to be perfectly periodic with the period T c \ = 27r in atomic units. The 
autocorrelation oscillates with the wave packet. The wave packet does not collapse, 
and there are no revivals. 

The autocorrelation for the particle in an infinite square well is shown in Fig. 5c. 
Here, T c \ ~ 0.04 a.u. and t rcv ~ 1.27 a.u. The initial periodicity of the wave packet 
is revealed by the periodicity in the peaks of |^4(t)| 2 near t = 0, with period T c \. The 
peaks decrease with time, corresponding to the collapse of the wave packet. Fractional 
revivals are evident at |t rev — 0.32 a.u., \t rew ~ 0.64 a.u., and |t rev — 0.95 a.u. The 
periodicity of |^4(t)| 2 at |t rev and |t rev is \T C \ and is evident in the graph in that the 




1 



-(p-po) 2 
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peaks are much closer together at these times. The autocorrelation |v4(t)| 2 = 1 at 
t rev , showing that the full revival is a perfect revival. The periodicity near t rev is equal 
to T c \ again. This cycle of periodicity, collapse, and fractional/full revivals repeats at 
multiples of t rev . There are no superrevivals. 

Figure 5d is the autocorrelation for the circular Rydberg wave packets considered 
in Sec. Ill, with the time given in nanoseconds. The initial periodicity T c \ ~ 0.26 nsec 
is difficult to resolve on the scale used in this graph. However, the full revival at t rev ~ 
21.0 nsec as well as some of the fractional revivals are evident. The large peaks at t m 
|t rev ~ 10.5 nsec indicate the formation of a single wave packet. This revival, however, 
is one half cycle out of phase with the classical motion. The full revival at t ~ t mv 
is evidently not a perfect revival since |^4(£)| 2 is less than 1. Fractional superrevivals 
appear at t ~ j^t sr ~ 105 nsec and t ~ j^t ST ~ 158 nsec with periodicities |t rev 
and \t revi respectively. The superrevival occurs at t « |t sr ~ 315 nsec, and the 
periodicity is \t tev . The height of the peaks at t ~ |t sr are greater than those at t rev 
indicating that the superrevival resembles the initial wave packet better than does 
the full revival at t m t rev . 

V. SUMMARY 

In this paper, a general treatment of the evolution and revival structure of local- 
ized quantum wave packets is presented. We consider wave packets formed from a 
superposition of states heavily peaked around a mean value n of the relevant quan- 
tum number, as would occur in an excitation process using a short-pulsed laser. The 
gaussian form for the weighting coefficients in Eq. (0) is used to simplify the calcula- 
tions in the examples, but the results do not depend on this form. Three time scales 
T ch t rev , and t sr controlling the evolution and revivals of wave packets are introduced. 
They are defined in Eq. (|5]) and are governed only by the dependence of the energy 
E n on the quantum numbers n. 

We analyze the revival structure of localized quantum wave packets for different 
quantum systems. Given the form of the energy spectrum E n for a given quantum 
system, it is possible to predict the motion and types of revivals that a wave packet 



14 



will undergo. We show that wave packets in certain classes of quantum systems 
exhibit particularly simple revival structures, with perfect classical periodicity or 
perfect full revivals. These results are illustrated with several examples taken from 
quantum systems that are commonly used in classes at the advanced undergraduate 
or first-year graduate level and that form the foundation for much of the analysis of 
wave packets in the literature. Plots are provided showing the wave packets and the 
autocorrelation functions for these examples. The revival structures analyzed here 
are applicable to a multitude of physical systems using wave packets in physics and 
chemistry. 
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Fig. 1: Unnormalized quantum wave packets for the simple harmonic oscillator with 
n — 15 and o = 1.5. The probability density is plotted as a function of x in 
atomic units at the times: (a) t — mod T d , (b) t = \T C \ mod T d , (c) t = \T C \ 
mod T ch (d) t = §T d mod T d . 

Fig. 2: Unnormalized quantum wave packets for a particle in an infinite square well 
with n — 15 and <r = 1.5. The probability density is plotted as a function of 
x in atomic units at the times: (a) t — mod t rcv , (b) £ = |T d mod t rev , (c) 
t = \T A mod t rcv , (d) t = \T cX mod t rcv , (e) t = T d mod t rcv . 

Fig. 3: Unnormalized quantum wave packets for the rigid rotator with I — 15 and 
a = 1.5. The probability density is plotted as a function of in radians at 
the times: (a) t — mod t rcv , (b) t = |T d mod t rcv , (c) t = T d mod t rcv , (d) 
t = \t rcv mod t rcv , (e) t = |t rcv mod t rcv , (f) t = it rcv mod t rev . 

Fig. 4: Unnormalized circular wave packets with n = 72 and <r = 2.5. The probability 
density is plotted as a function of x and y in atomic units at the times: (a) t — 0, 
(b) t = |T d , (c) t = T ch (d) t = \t Tew , (e) t = t rev , (f) i = |t sr . 

Fig. 5: The absolute square of the autocorrelation function as a function of time for 
four different quantum systems: (a) the free particle, (b) the simple harmonic 
oscillator, (c) the infinite square well, (d) Rydberg wave packets. 
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